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ABSTRACT 



a. 

^ , Context. Mercury-manganese (HgMn) stars have been considered as non-magnetic and non- variable chemically peculiar (CP) stars 

for a long time. However, recent discoveries of the variability in spectral line profiles have suggested an inhomogeneous surface 
' distribution of chemical elements in some HgMn stars. From the studies of other CP stars it is known that magnetic field plays a key 

I role in the formation of surface spots. All attempts to find magnetic fields in HgMn stars have yielded negative results. 

Aims. In this study, we investigate the possible presence of a magnetic field in (p Phe (HD 1 1753) and reconstruct surface distribution 
of chemical elements that show variability in spectral lines. We also test a hypothesis that a magnetic field is concentrated in chemical 
spots and look into the possibility that some chemical elements are stratified with depth in the stellar atmosphere. 
Methods. Our analysis is based on high-quality spectropolarimetric time-series observations, covering a full rotational period of 
the star. Spectra were obtained with the HARPSpol at the ESO 3.6-m telescope. To increase the sensitivity of the magnetic field 
^ ^ i search, we employed the least-squares deconvolution (LSD) technique. Using Doppler imaging code INVERSIO, we reconstructed 

^ f surface chemical distributions by utilising information from multiple spectral lines. The vertical stratification of chemical elements 

was calculated with the DDAFit program. 

Results. Combining information from all suitable spectral lines, we set an upper limit of 4 G on the mean longitudinal magnetic field. 
For chemical spots, an upper limit on the longitudinal field varies between 8 and 15 G. We confirmed the variability of Y, Sr, and Ti 



o 

^ ■ and detected variability in Cr lines. Stratification analysis showed that Y and Ti are not concentrated in the uppermost atmospheric 

^ ' layers. 

Conclusions. Our spectropolarimetric observations rule out the presence of a strong, globally-organised magnetic field in (p Phe. 

. This implies an alternative mechanism of spot formation, which could be related to a non-equilibrium atomic diff'usion. However, the 

^ ■ typical time scales of the variation in stratification predicted by the recent time-dependent diff'usion models exceed significantly the 

I spot evolution time-scale reported for if Phe. 



o 
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1. Introduction global magnetic field sufficient to cause the chemical spot for- 

' ! mation. 

Mercury-manganese (HgMn) stars is one of the sub-classes of ^ ^- ^ ^ • ^ • ^ ^ ^ ^ u at 

^-H 1 , . 11 1. /^T^x TT Tv^r 1 Detection or spots mcreased the mterest towards HgMn 

early-type chemically peculiar (CP) stars. HgMn stars show an ^ i^- • r j u a/t 

! 11 i^TT tT/t o 1 1 11 1 stars, resulting m a discovery or seven more spotted HgMn 

> overabundance of Hg,Mn,Y,Sr and other, mostly heavy, chem- h r ^SS and HR 8723 (iKo chukhov et alj 1^. 

^ ical elements. These stars are often found in close binaries and faubrig et alJ 1200^ iFolsomlTa l. 2010), HD 1 1753, 

^ . belong to the temperature range of Teff = 9500-16 000 K, cor- hD 53244, HD 22 1507 dBriquet et alJ r2010). and HD 32964 

respo nding to the spectral classes from AO to B5 (IDworetskvl jrn r-j — — j-i \ r TTIZ ^- tt 

i .Q^^ L ^ ^ ^ ^ (IMakaganiuk et al.l 1201 lal) . At the same time, all systematic 

' ^' studie s of the magnetic field in HgMn stars reported its ab- 

Earlier it was thought that HgMn stars show no rotational ^ence ( Shorlin et ai;2002':'Wade et ai:2006':'Fbli m^r^rd^ 

spectral variability, meaning that the atmospheres of these stars Auriere et al. 2010; Makasaniuk et al. 201 lb). In the best cases, 

are horizontally homogeneous. However, the variability discov- ^n upper limit set on the strength of the mean longitudinal mag- 

ered in a spectral line of Hg in a And (Adelman et al. 2002) netic fi eld is a few Gauss (Auriere et ZI[2QTqI:Im^ saniuk et al. 

has changed this picture. I Adelman et al.l (120021) show that the |2Qiib|) . 

observed variability is related to presence of chemical spots. In .^tta inco tttt^ oooo ttt^ cco\ • 

,1 r A . .1 r ^ ' 11 11 The source w Phe (HD 11753, HIP 8882, HR 558) is a 

the case ol Ap stars, the presence or sp ots is usu ally caused by , , ^ ^. T . i , c^^\ j ^r^r^r^r^ 

^1 ^. 2 ,1 rj.. /i 1 ,^.1 l~nJRKR^\ ' A slowly-rotating, bright (V = 5.11) an d cool {Lq,^- 10700 K , 

the magnetic field The study by Wade eF^ (|2006) aimed at he [s^iith & pZretsklTQQl HgMn star, bworetskv et al.l dlOSl 

search for magnetic field m a And proved an absence of the TTirTT^ — — \ v r t^i^ ^7^ . , . 

^ ^ suggested that the radial velocity ol (f Phe varies with a period or 

more than 30 days, which me ans that this star is a s pectro scopic 



more tnan 30 days, wnicn me ans tnat tnis star is a s pectro scopic 
binary with a single spectrum. iLeone & Catanzarol (119991) found 



* Based on observations collected at the European Southern 
Observatory, Chile (ESO programme 084.D-0338). the orbital period of 41 .489 ± 0.019 d for (f Phe. 
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Emp loying the lUE spectra a n d the spectrum synthesis 
method, "Smith & Dworetsy (Il993h : ISmithI (Il993l 11994 Il997h 
determined a bundances of Cr, Mn , Fe, Mg, Al, Si, Zn, C u, Co, 
Ni, and Hg. IJomaron et alJ (1 19991) and iDolk et alJ (l2003h used 
optical spectra to determine abundances of Mn and Hg, respec- 
tiv ely. 

ISmithI (Il996l) found a stratification of gallium, based on the 
analysis of the UV lines of this element. He concluded that Ga 
shows an increase of concentration, starting from logrsooo = 
+0. 3 towards the upper layers of the stellar atmosphere. 

iBriquet et al.l (l2QlQh discovered chemical spots on ^ Phe. 
Using a large number of observations, the authors found vari- 
ability in the spectral lines of Ti, Sr, and Y with a 9?54 period. 
They reconstructed surface maps for the two sets of spectra of 
this star, separated by 65 days. Based on the diflTerences in their 
maps, the authors suggest an evolution of spots, si milar to the 
phenomenon discovered by lKochukhov et af] (l2007h in a And. 

There were no previous magnetic field studies of ip Phe. 
This star was included in the sample of HgMn stars for 
which we perfornied a m agnetic field survey with HARPSpol 
(iMakaganiuk et al.l l2Qllb ). The spectropolarimetric observa- 
tions of (f Phe cover its full rotational period, enabling us to 
measure the magnetic field at each rotational phase. Based on 
these high-quality data we also reconstructed surface maps of 
variable chemical elements and tested the presence of extreme 
vertical stratification for some of them. 

In Sect. [2] we describe the observations and data reduction. 
The least- squares deconvolution (LSD) and the measurements of 
the longitudinal magnetic field are presented in Sect.O Section|4] 
describes investigation of the line profile variability. Surface 
maps of chemical elements derived with the Doppler imaging 
technique are discussed in Sect. [S] The stratification analysis is 
presented in Sect. [6l Section [7] summarises our results and dis- 
cusses them in the context of previous studies of (p Phe and re- 
cent theoretical diflTusion calculations. 

2. Observations and data reduction 

The star (p Phe was observed i n January 2010 , using a newly- 
built HARPS polarimeter (Snik et al.1 120111: iPiskunov et aP 
1201 ll) at the 3.6-m ESO telescope in La Silla, Chile. All obser- 
vations were done with the circular polarimeter. Data collected 
during the first two nights were part of the engineering run de- 
voted to the commissioning of the instrument. Orientation of the 
quarter- wave retarder plates was not calibrated, therefore no use- 
ful Stokes y spectra could be obtained. During the ten other 
nights, we obtained complete spectropolarimetric observations 
of (p Phe. A typical signal-to-noise {S IN) ratio of our spectra is 
300-400 per CCD column dX A ^ 5200 A and the resolving 
power is A|^A= 115000. 

Spectra were recorded by a mosaic of two 2K x 4K CCDs. 
One CCD records the blue part of a spectrum, which contains 
90 spectral orders in the spectropolarimetric mode, while the 
red part of a spectrum is recorded by the other CCD. The lat- 
ter contains 52 spectral orders. The overall wavelength coverage 
is 3780-6920 A, with a gap between 5259 and 5337 A. For the 
calibration procedure we used a standard set of frames, which 
consisted of ten biases, ten flat fields exposed in circular polar- 
isation mode, and one ThAr frame. Such a set was obtained at 
the beginning and end of each observing night. The polarimet- 
ric observations consisted of four sub-exposures, corresponding 
to 45°, 135°, 225° or 315° position of the quarter-wave plate. 
The exposure times adopted for individual sub-exposures varied 
from 270 to 300 seconds, depending on the weather conditions. 



Table 1. The log of the HARPSpol observations of (p Phe. 



HJD-24 X 10^ 


Phase 


Stokes 


^exp (s) 


SIN 


55197.06789 


0.000 


I 


4x48 


190 


55198.04517 


0.102 


I 


4x47 


183 


55198.05107 


0.103 


I 


4x200 


507 


55198.06359 


0.105 


I 


4x300 


388 


55199.04987 


0.208 


I 


4x200 


352 


55199.07179 


0.210 


I 


4x41 


248 


55199.09245 


0.212 


I 


4x200 


365 


55200.05063 


0.313 


IV 


4x200 


299 


55201.02209 


0.415 


IV 


4x200 


395 


55202.03895 


0.521 


IV 


4x200 


433 


55203.02838 


0.625 


IV 


4x300 


348 


55204.01590 


0.728 


IV 


4x280 


430 


55205.04475 


0.836 


IV 


4x280 


353 


55206.02121 


0.939 


IV 


4x270 


406 


55207.01430 


0.043 


IV 


4x270 


476 


55209.00910 


0.252 


IV 


4x270 


379 


55210.03112 


0.359 


IV 


4x270 


564 



The rotational period of Phe, 9j54 (iBriquet et al.ll2010l) . 



was recently revised to 9?531 (iKorhonen et al.l 



120111) . This im- 



provement corresponds to a negligible change of the rotational 
phases of o ur observa tions. Therefore, we adopted the period 
published by lBriqueret al.l (|2010) throughout our analysis. In the 
Table [T] we list the heliocentric Julian dates of the observations, 
rotational phase, type of the observations, exposure times, and 
S IN. The rotational phases were computed relative to the first 
night of observations, assuming the rotational period of 9?54 and 
the zero phase at HJD = 2455197.06789. 

To process our data w e employed the REDUCE package 
dPiskunov & Valentill2002h . which performs the standard reduc- 
tion steps, such as averaging of bias and flat field frames and 
subtracting the master bias from the master flat field and sci- 
ence frames. Each science frame was corrected for the pixel- 
to-pixel sensitivity variations using a normalised flat field. All 
stellar spectra were corrected for the scattered light before the 
optimal extraction. 

The wavelength calibration was performed based on one 
ThAr spectrum for each night, yielding an accuracy of 18- 
21 m s"^ of the derived wavelengths. Using a single comparison 
spectrum per night is sufl&cient for our study, thanks to the in- 
trinsic stability of the HARPS spectrometer (Mavo r et al.ll2003l) . 
The normalisation of spectra was done first by correcting the 
blaze shape and then dividing by a smooth function that was fit- 
ted to an upper envelope of the merged spectrum. 

To minimise the possible eflfects of spurious polarisation, we 
employed the ratio method (Bagnulo et al .1200 9) to combine po- 
larimetric sub-exposures and derive continuum-normalised cir- 
cular polarisation (V//c) spectra. For brevity, we refer to these 
data as "Stokes V spectra". In addition to the polarisation spec- 
tra, the ratio method produces a diagnostic spectrum, which is 
called the null (N) spectrum. To derive the null spectrum, the 
four sub-exposures are combined in a such way that the stellar 
polarisation is destroyed, enabling us to diagnose possible spu- 
rious polarisation. 

3. Magnetic analysis 

3.1. LSD technique 

A search for weak magnetic fields requires polarimetric data of a 
very high S IN. Achieving this directly is technically challenging 
and costly in terms of the observational time. On the other hand. 
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the circular polarisation profiles are very similar in shape for the 
majority of spectral lines, especially in weak fields. Therefore, 
an alternative approach to increasing polarimetric sensitivity 
while maintaining reasonable exposure times is to use a multi- 
line diagnos tic technique, such as the least- squares de convolu- 
tion (LSD) ( Donati et al.lll997HKochukhov et al.ll201Qh . 

This technique assumes that all spectral lines are identical 
in shape and can thus be described by a scaled mean profile. If 
we assign a certain set of weights to individual spectral lines, 
then the mean profile multiplied by these weights should re- 
produce an observed spectrum. Utilising observations and a set 
of weights, LSD code computes a mean profile, combining in- 
formation from many spectral lines. As a result, the S IN in- 
creases by a factor of up to 30-40, depending on the number 
of lines in a stellar spectrum. LSD has pro ven to be a power- 
ful tool for detecting weak magnetic fields ( Auriere et al.ll2009l: 
iLignieres et al.l 12009V For HgMn stars, it allowed for an upper 
limit to be established f or the mean longitudinal magnetic field 
at the level of < 10 G (I Auriere et al.ll2010l: iMakaganiuk et al.1 
I2011bh in the best cases. 

3.2. LSD mask and stellar parameters 

LSD technique requires a list of spectral line positions, their cen- 
tral intensities, and magnetic sensitivity, collectively called a line 
mask. To compile such a mask, it is important to adopt correct 
stellar parameters since the intensities and the total number of 
spectral line s extracte d from t he Vienna Atomic Line Database 
(VALD, Pis kunov et al. 1995; Kupka et al. 1999) will depend on 
the eflTective temperature and chemical composition of a star. 

We estimated Teff and logg from the Stromgren photome- 
try ( Ha uck & Mermimodl[T998h with the he lp of the TempLogG 
code ( Kaiser"2006'), using the calibration by Moon & D woretskyI 
([l985). The resulting T^q = 10500 K and log ^ = 3 .8 agree with 
the parameters determined by Smith & Dwo retskyI ([1993) (Jeff = 
10700 K, logg = 3.8) and .Dolk et al. (2003^ (T^^ = 10612 K, 
logg = 3.79). The scatter of the eflTective temperature deter- 
minations suggests Teff uncertainty around 200 K. According to 
|Moon & Dworetsky (1985), the uncertainty of the eflfective tem- 
perature and surface gravity established using their calibration 
is +260 K and ±0.1 dex, respectively. The mod el atmosphere 
of (fi P he was computed with the LLmodels code (IShulyak et al.l 
l2004j) for r.ff = 10500 K and logg = 3.8. 

The abundances of many chemical elements in (f Phe are un- 
known, although some other elements have been studied previ- 
ously (see Sect.[T]). Using these individual abundance and adopt- 
ing average values typical of HgMn stars for elements not anal- 
ysed before, we compiled a preliminary lin e list and com puted 
synthetic spectrum with the SYNTH3 code (lKochukhovl [2007). 
Chemical abundances were refined by comparing this calcula- 
tion to the observed spectrum obtained at HJD = 55210.03112, 
which has the highest S /N. Zero microturbulent velocity was 
assumed both for the VALD extraction and for the spectrum syn- 
thesis. 

We found good agreement between our observations and 
the synthetic spectrum for the spectral lines of chemical ele- 
ments studied previously. However, since we adopted a some- 
what cooler model atmosphere compared to lSmith & DworetskyI 
(1993), we had to reduce the abundances of Fe, Cr, Mg, Mn, and 
Ga by up to ^ 0.4 dex. No changes were required for Si and Co. 
Relatively large corrections were introduced for Sc, Y, Pt, and Zr, 
because spectral lines of some of these chemical elements were 
excessively strong or too weak, mismatching the observations. 
Since the precise quantitative abundance analysis is beyond the 



Table 2. Abundances of chemical elements in (f Phe. 



Element 


Initial abundance 


Final abundance 


He 


-2.07''' 


-1.65 


N 


-4.26^ 


-4.46 





-3.58''' 


-3.45 


Na 


-5.87''' 


-5.35 


Mg 


-4.64^ 


-5.00 


Al 


-7.20^ 


-7.20 


Si 


-4.68^ 


-4.70 


P 


-5.68^ 


-6.00 


S 


-5.40^ 


-5.20 


Ca 


-5.53"^ 


-5.73 


Sc 


-8.49^ 


-10.75 


Ti 


-6.64^ 


-6.50 


V 


-8.04^ 


-9.54 


Cr 


-5.80^ 


-6.00 


Mn 


-5.80^ 


-5.90 


Fe 


-4.35^ 


-4.45 


Co 


-9.50^ 


-9.50 


Ni 


-6.90^ 


-6.90 


r^n 
L.U 


— 1 .jL 


— 1 .dL 


Zn 


-lO.OO^ 


-10.00 


Ga 


-7.10^ 


-7.30 


Sr 


-7.12"^ 


-7.35 


Y 


-7.83"^ 


-6.95 


Zr 


-7.45^ 


-8.50 


Pt 


-7.66^ 


-6.75 


Hg 


-6.88^ 


-6.88 



Notes. References to the initial abundance are marked as follows: 
average abundance estimate for HgMn star, Smith & Dworetsk\i 
(1993) , Smith (1993) ^^^ISmithl (^994), ^''^ .Smith (.1996) . .Smith! 
(1199% . ^^^ iJomaron et all (119991) . 



scope of our study, we adjusted abundances of individual chemi- 
cal elements visually, without deriving formal error bars. This is 
sufficient for compiling an LSD line mask. The abundances are 
given in the Table [2] as follows: chemical element and the initial 
and final abundances. The revised abundances were used to ob- 
tain a new line list from VALD, which provided 3827 spectral 
lines in the 3781-6915 A range. 

To apply LSD, we constructed a line mask containing central 
wavelengths of spectral lines and a set of weights for Stokes / 
and y. Central depths returned by VALD were used as weights 
for the intensity spectrum. For the circular polarisation spectrum, 
the weights were calculated as wy = dAg/Ao, where d is the 
central depth of a spectral line, A the wavelength, and g the cor- 
responding eff'ective Lande factor. The normalisation parameter 
Aq = 4696 A was chosen to be close to the averaged wavelength 
of all spectral lines in the mask. 

To get optimal results from the LSD technique, it is essential 
not to include very weak spectral lines, because an extra noise 
might be added to the LSD profiles. To exclude weak lines, we 
apply a cutoff' criterion, which sets a threshold for the line inten- 
sity. All spectral lines weaker than the cutoff* level are excluded 
from the LSD mask. For (p Phe we se t the cutoff criterion to 
0.05, which means that our LSD code ( Kochukh ov et al 
will use only those lines that are deeper than 5% relative to the 
continuum level in the unbroadened spectrum. This yielded 962 
spectral lines for computing LSD profiles. 
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Table 3. Magnetic field measurements of if Phe. 



0.08 



, 0.04 



0.00 
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Fig. 1. LSD profiles of (f Phe for Stokes / (left panel) and 
Stokes V (right panel) spectra. For display purposes we shifted 
profiles vertically. The rotational phase for each profile is de- 
noted on the left. 



3.3. Analysis of LSD Stokes V profiles 

Taking Vrad = 14 kms"^ of (p Phe into account, we recon- 
structed LSD profiles for the velocity range between -36 and 
64 kms"^ which makes the profile extension symmetric with 
respect to the line centre. An average pixel of the HARPS CCDs 
has a resolution of 0.8 km s"^ We adopted this value as a step for 
LSD profiles. The S /N reported for the final LSD profiles cor- 
respond to this velocity bin. The average S /N gain is nine. The 
LSD Stokes / and V profiles for ten rotational phases are shown 
in Fig.[T] Based on visual examination of the LSD V profiles, we 
conclude that no polarisation signal is present. 

We estimated the mean longitudinal magnetic field,(5z), 
from the first mome nt of the LSD Stokes V profile 
(iKochukhov et al.|[201 (j). setting the mean Lande factor and cen- 
tral depth to unity, in accordance with the adopted normalisa- 
tion of LSD weights. All profiles were integrated from -4 to 
34 km s"^ These limits were chosen symmetrically with respect 
to the centre of the LSD / profile and were utilised for (5z) mea- 
surements in the LSD V and LSD N profiles. The use of the 
null profile is an important diagnostics of <5z) measurements in 
LSD V. Should a spurious polarisation alter Stokes V spectrum, 
it will be indicated by the (5z)(N) measurements significantly 
diflTerent from zero. 

Our analysis of the (5z)(N) measurements showed the ab- 
sence of a spurious polarisation signal. At the same time, none 
of (5z)(V) measurements indicates there is a magnetic field in 
(f Phe. The average error of (5z) is 3.8 G and the minimum is 
2.8 G. 

In addition to the longitudinal field measurements, we per- 
formed a false alarm probability (FAP) analysis. This method, 
based on statistics, diagnoses a statistical significance of the 
deviation of the LSD profile from zero. The FAP analysis is espe- 







LSDV 




LSDA^ 




Phase 


S/N(LSD) 


(5z> (G) 


FAP 


(5z> (G) 


FAP 


0.043 


2532 


-9.7 + 5.1 


1.00 


-9.0 + 5.1 


1.00 


0.252 


4069 


-3.5 + 3.2 


0.66 


-5.4 + 3.2 


0.89 


0.313 


3130 


6.6 + 4.0 


0.25 


4.4 + 4.0 


0.49 


0.359 


4357 


1.3 + 2.9 


0.60 


3.7 + 2.9 


0.00 


0.415 


4417 


1.9 + 2.8 


0.16 


-3.7 + 2.8 


0.23 


0.521 


3053 


-6.0 + 4.1 


0.78 


-1.4 + 4.1 


0.68 


0.625 


3404 


6.8 + 3.8 


0.80 


-1.6 + 3.8 


0.47 


0.728 


3311 


-1.3 + 4.0 


0.36 


1.8 + 3.9 


0.66 


0.836 


3954 


3.3 + 3.3 


0.78 


2.9 + 3.3 


0.32 


0.939 


4482 


-0.7 + 2.9 


0.97 


-3.8 + 2.9 


0.26 



o 

CD 




0.2 0.4 0.6 0.8 
rotational phase 

Fig. 2. Measurements of the longitudinal magnetic field as a 
function of rotational phase. 



cially sensitive to complex magnetic field configurations, which 
may yield zero longitudinal field. We interpreted FAP results fol- 
lowing Donati et aL| (119971) . Any FAP larger than 10"^ is consid- 
ered as no detection of the magnetic field. Our FAP values are 
always greater than 0.1, from which we conclude that (f Phe sig- 
nificant tangled magnetic fields. 

We summarise the longitudinal magnetic field and FAP mea- 
surements in Table[3l The first column gives the rotational phase, 
the second column gives S/N in the LSD Stokes V profiles, 
while the last four provide longitudinal magnetic field measure- 
ments for LSD V and N with the corresponding FAP numbers. 

To check for a systematic rotational modulation of (5z), we 
plotted our measurements as a function of the rotational phase 
(Fig. O. We conclude that there is no evidence of periodicity. 
All measurements are consistent with zero longitudinal magnetic 
field. 

It has been suggested by 'Hubrig et al.' ('2010') that magnetic 
fields in HgMn stars may be concentrated in spots of chem- 
ical elements, implying that the polarisation signal should be 
stronger in the variable spectral lines. To constrain the longitudi- 
nal magnetic field in chemical spots, we computed mean profiles 
for Y, Ti, and Cr using the multi-profile version of our LSD code 
(see Sect. 2.4.1 in Kochukhov et al. 2010). From these LSD pro- 
files, we measured (5z), finding no magnetic field with a typical 
error bar of 8 G for Ti, 13 G for Cr, and 15 G for Y lines. These 
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0.4 0.6 0.8 

rotational phase 

Fig. 3. Measurements of the longitudinal magnetic field as a 
function of rotational phase, based on spectral lines of Y (cir- 
cles), Ti (triangles), and Cr (squares). 

measurements are plotted as a function of rotational phase in 
Fig.H 

Summarising magnetic field analysis results for the inhomo- 
geneously distributed chemical elements, we conclude that there 
is no evidence of any polarisation signal in their lines. 




5.51966-104 5.52003-104 5.52041-104 5.52078-104 5.52115*104 

HJD 




5.5198-104 5.5209-104 



5.5220-104 
HJD 



5.5231-104 5.5242-104 



Fig. 6. Upper panel: measured radial velocity with respective 
error bars. Lower panel: the orbital solution (solid line) based 
on the orbital elements from Leone & Catanzaro (1999) and our 
measurements (circles) as a function of reduced Julian date. 



4. Line profile variability 

For the line profile variability search, we utilised all available 
intensity spectra from our observing run. We examined the spec- 
tral range from from 3830 to 6830 A, ignoring all stellar spectral 
lines blended by telluric absorption. In total, we found 83 vari- 
able spectral lines. 

The strongest variability is seen in the lines of Y and Sr. 
Distortion of the shape of these lines is clearly evident in indi- 
vidual spectra. A somewhat weaker variability is observed for Ti. 
We see a diflTerent variability in spectral lines of various chemical 
elements, but a similar pattern for the lines of the same element. 
This type of variability is typical of the spotted distribution of 
chemical elements. No definite variability was found in individ- 
ual spectral lines belonging to elements other than Y, Sr, and Ti. 
Compared to other CP stars, the HgMn type often shows a rather 
weak variability, which is barely detected based on the visual ex- 
amination of individual spectral lines. Thus, a numerical assess- 
ment and multi-line analysis are necessary for investigating this 
weak variability. 

We computed a standard deviation for all spectral lines. 
Examination of these data confirmed the variability in Ti, Sr, 
and Y, but also indicated a marginal variability in Cr. We show 
17 phases for selected spectral lines of Ti, Cr, Sr, and Y in Fig.lH 
The uppermost curve represents the standard deviation, which 
shows a distinctive profile at the position of variable spectral 
lines. The amplitude of the standard deviation profiles charac- 
terises the strength of variability. 

We applied the LSD technique to the spectral lines of Cr to 
verify their variability. Results of this analysis are summarised 
in Fig. [5] For the reference, the LSD profiles for the non- 
variable chemical element, Fe, are shown in the left-hand panel. 
The other two panels present LSD profiles for Cr and Y. We 
computed an average LSD profile for each chemical element 



and plotted it with dashed lines to highlight the relative pro- 
file changes. Indeed, the very weak variability suggested by the 
standard deviation is seen in the LSD profiles of Cr. The com- 
parison of Fe and Cr profiles allows us to conclude that these 
weak variations are not introduced by possible instrumental or 
data-reduction artefacts. 

We m entioned in Sect . [T] that 'Dwor etskv et al.1 (Il982h and 
iLeone & Catanzarol (Il999h have suggested that the radial ve- 
locity of if Phe changes due to the orbital motion. To verify 
their results, we used the centre-of-gravity of the LSD / pro- 
files of Fe to determine the radial velocities for each observa- 
tion. The resulting mean value of the radial velocity of if Phe is 
14. 58+0. 02km s"^ Individual measurements of radial velocity, 
together with the respective error bars, are shown in the upper 
panel of Fig. [6l T he radial velocity va riation predicted by the 
orbital solution by llLeone & Catanzarol ([1999) is compared with 
our measurements in the lower panel of Fig. O The errors of 
individual measurements are close to the symbol size. 

We conclude that the radial velocity of the star does not 
change with the suggested 41^489 period. If this period were 
true, we would see a change in radial velocity ^^10 kms"^ dur- 
ing our o bserving run Moreover, earlier observations of if Phe 
by Br iquet et al.l (1201 Oh with a better time coverage also indicate 
no change in the radial velocity of this star, yielding the mean 
radial velocity, 14.24+0.01 kms"\ nearly identical to ours. 

We use Y lines, which show the strongest variability com- 
pared to other chemical elements, to verify the rotational period 
of 9?54 suggested for (p Phe by B riquet et al. (2010). The radial 
velocities derived from the LSD / profiles of Y (see Fig. O are 
shown as a function of rotational phase in Fig. [71 Evidently, our 
data agrees with the previously determined rotational period. A 
short time span of our observations does not allow deriving any 
useful independent period estimate. 
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Fig. 4. The variability of spectral line profiles in if Phe. The uppermost curve is a standard deviation computed for all spectra below. 
For displaying purposes, spectra are shifted vertically. 
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Fig. 7. Radial velocities determined from the LSD / profiles of 
Y as a function of rotational phase. 



5. Doppler imaging 

A quantitative interpretation of the observed variations in spec- 
tral lines requires using the Doppler imaging (DI) technique, as 
well as knowledge of physical parameters of a star. One of the 
key parameters is the projected rotational velocity. To determine 
Ve sin /, we selected 20 sufl&ciently strong unblended spectral 
lines of iron in the spectrum with the highest S /N (see Table [T]). 
We fitted the model spectrum with corrected abundances to these 



observations, using BinMag^. BinMag is an IDL-based widget 
tool, which allows visual comparison of synthetic spectra with 
observations, taking radial velocity, instrumental resolution, and 
rotational and macroturbulent broadening into account. To de- 
termine the rotational velocity of (f Phe, the code was allowed to 
vary Ve sin / for each selected Fe line. The resulting mean value is 
Ve sin / =13.62 km s"^ with a standard deviation of 0.22 km s"^ . 
We have verified that variation by one standard deviation yields 
spectrum changes above the noise level, therefore this error esti- 
mate is meaningful. 

Another parameter required by DI is the inclination angle / 
between the stellar rotation axis and the line of sight of an ob- 
server. It can be determined as 



smi 



P Ve sin / 
50.6137?' 



(1) 



where P is the stellar rotational period in days, Ve sin / the pro- 
jected rotational velocity, and R the radius of a star in solar units. 
The last is given by the following expression: 



(2) 



The stellar luminosity can be inferred from the absolute magni- 
tude: 



L\ My + BC- Mf 



bol 



2.5 



(3) 



Here My is an absolute visual magnitude of the star, BC is the 
bolometric correction, and .73. The bolometric correc- 

tion for (f Phe was estimated following Floweij (Il996h . The ab- 
solute visual magnitude was determined as My = V + 5 + log tt. 



http://www.astro.uu.se/~oleg/soft.html 
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Fig. 5. Constant LSD profiles of Fe (left panel), extremely weakly variable Cr profiles (middle panel), and LSD profiles of Y 
(right panel) as an example of the strongest variability. Dashed line indicates the mean profile. For displaying purposes, profiles 
corresponding to diflferent rotation phases are shifted vertically. 



Table 4. Physical properties of ip Phe. 



Parameter 


Value 


/ 


65.7° +7.r 


M, 


0.243 + 0.076 


R/Rq 


2.817 + 0.157 


log L/Lq 


1.938 + 0.035 


logg 


3.8 + 0.1 


Teff (K) 


10500 + 200 


Ve sin / (kms~0 


13.62 + 0.22 



The t rigonometric p arallax, tt = 10.63+0.37 mas, was adopted 
from Ivan Leeuwenf 12007). The eflTective temperature is taken 



from Sect. 13.21 From all these parameters we derive the incli- 
nation angle / = 65.7° + 7.1°. All parameters of (p Phe and their 
respective errors are summarised in Tabled 

To interpret the variability of line profiles, we selected 
29 spectral lines of Y, Sr, Ti, and Cr. Surface maps were 
reconstructed using the Dop pler imaging code INVERSIO 
dPiskunov & Kochukhovll2002b . Unlike the previous attempt to 



study chemical spots in (f Phe (iBriquet et al.ll201Ql) . our anal- 
ysis is based on simultaneous modelling of multiple spectral 
lines, which significantly improves reliability of the DI maps. 
The maps of Y were reconstructed from nine spectral lines, Sr 
from three, Ti from nine, and Cr from eight spectral lines. 

All inv ersions were constrained using the Tikhonov regular- 
isation (see'Piskunov & Kochukhovl2002h . The value of the reg- 
ularisation parameter was chosen in the usual way, by adopting 
the highest regularisation which provided a fit to the observa- 
tions compatible with S /N of the data. 

Analysing the fit of individual profiles, we found synthetic 
spectrum deviating systematically for a couple of lines of each 
element, while the rest of lines were fitted rather well. To min- 
imise the influence of the errors in atomic data, we allowed the 
code to correct the oscillator strengths of these lines in the course 
of inversion. In Table \5\ we provide the information about ion, 
wavelength, excitation potential, initial oscillator strength, and 
corrected oscillator strength of the spectral lines used in DI in- 
version. 
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Table 5. Spectral lines used for the Doppler imaging of ip Phe. 



Ion 


Wavelength (A) 


Excit (eV) 


Initial \oggf 


A logg/ 


Yii 


4204.695 


0.000 


-IA61 


-0.226 


Yii 


4398.013 


0.130 


-0.895 




Yii 


4465.268 


4.133 


0.222 


-0.195 


Yii 


4900.120 


1.033 


0.103 




Yii 


5087.416 


1.084 


-0.077 




Yii 


5119.112 


0.992 


-1.249 




Yii 


5497.408 


1.748 


-0.276 




Yii 


5509.895 


0.992 


-0.948 


0.076 


Yii 


5662.925 


1.944 


0.384 


0.167 


Srii 


4077.709 


0.000 


0.160 




Srii 


4215.519 


0.000 


-0.140 


0.057 


Srii 


4305.443 


3.040 


-0.130 


-0.355 


Tin 


4163.644 


2.590 


-0.130 




Tin 


4386.844 


2.598 


-0.960 




Tin 


4395.839 


1.243 


-1.930 




Tin 


4399.765 


1.237 


-1.190 




Tin 


4501.270 


1.116 


-0.770 


0.203 


Tin 


4563.757 


1.221 


-0.690 




Tin 


4779.985 


2.048 


-1.260 




Tin 


4911.193 


3.124 


-0.610 




Tin 


5129.152 


1.892 


-1.240 


-0.035 


Cm 


4037.972 


6.487 


-0.670 


-0.238 


Cm 


4051.930 


3.104 


-2.330 


-0.155 


Cm 


4145.781 


5.319 


-1.100 


-0.114 


Cm 


4275.567 


3.858 


-1.730 




Cm 


4554.988 


4.071 


-1.490 


0.052 


Cm 


4565.739 


4.042 


-1.980 




Cm 


4588.199 


4.071 


-0.840 


0.198 


Cm 


4616.629 


4.072 


-1.570 





We present spherical maps of the distribution of chemical el- 
ements in Fig. [HI The fit of the observed profiles of Y, Sr, Ti, 
and Cr is shown in Figs.[9l[T0l[TTl and[T2l respectively. For Srii 
line at 4077.71 A, we found a much stronger contribution of the 
blending spectral line of Cm at 4077.51 A in the synthetic spec- 
trum than appears in observations. The correction of the oscilla- 
tor strength for this particular Cr line by one order of magnitude 
allowed a better fit of the observed profile, reducing the standard 
deviation for Sr from 1 .06% to 0.65%. For the other chemical el- 
ements, the standard deviation of the final profile fit is 0.52% for 
Y, 0.45% for Ti, and 0.51% for Cr. These results are consistent 
with the quality of our data. 

From the DI surface maps we found that all chemical ele- 
ments concentrate in the same spot, observed at phases 0.2-0.4. 
A secondary chemical spot exists for Ti and Cr and is visible at 
phase range 0.8-1.0. However, the abundance gradients relative 
to the mean stellar composition are diflTerent for the four investi- 
gated elements. The largest gradient (1.3 dex) is found for Y. Sr 
shows less of a contrast (0.6 dex), which is comparable to that of 
Ti (0.5 dex). Finally, a very low contrast (0.15 dex) is observed 
for Cr. 

There is a significant variation in the latitudinal distribution 
of chemical elements. In particular, the relative abundances are 
diflTerent at the visible rotational pole. For Cr we observe an un- 
derabundance at the pole, whereas for Ti, Sr, and Y the over- 
abundance spots extend to high latitudes in a very similar way. 
The Cr underabundance zone occupies the whole polar region, 
while those of Ti, Sr, and Y occupy roughly half of it. 

One can also note an interesting variation in the overall sur- 
face area covered by the spots of diflTerent elements. A major 
part of the stellar surface shows a relative Ti overabundance. A 
somewhat smaller area is covered by Y spots, and even smaller 



occupied by Cr inhomogeneities. Sr is concentrated mainly in a 
single spot. 

6. Stratification of Y and Ti 

It is known that chemically peculiar stars may possess an in- 
homogeneous distribution (stratification) of some chemical el- 
emen ts with depth in stellar atmosphere (iRyabchikova et al.l 
[2QQ2h . In this case, computing a synthetic spectrum under the 
assumption of a homogeneous atmosphere will lead to discrep- 
ancies between the observed and synthetic profiles of spectral 
lines formed at diflTerent atmospheric levels. However, such a 
discrepancy must be interpreted with caution due to a possible 
confusion with an inadequate quality of the adopted line transi- 
tion data. 

A complete stratification analysis of (f Phe is beyond the 
scope of our paper. Instead, we aim at testing a very specific 
hypothesis that inhomogeneously distributed chemical elements 
are concentrated in the upper layers of the stellar atmosphere, 
where they might be affected by an extremely weak magnetic 
field. This sc enario is inspired by t he model of Hg stratification 
suggested by lMichaud et al.l (Il974l) and by the recent di scussion 
of the role of weak magnetic fields ("Alecian e t al .11201 ll) . We se- 
lected Y and Ti for the stratification analysis because their spots 
show a high relative abundance contrast, and these elements have 
a sufl&cient number of diverse spectral lines. The spectra for 
stratification analysis were obtained by averaging the line pro- 
files utilised for DI. 

We have carefully examined atomic data available for the 
selected Yii lines. Our primary source of transition probabili- 
ties is the revised Y ii list kindly calculated by R.Kurucfl These 
data are compatible with the oscillator strengths extracted from 
VALD, but include more lines. For seven out of nin e Y ii lines al - 
ternative log gf values are available from Biemont etal.1 (l 201ll) . 
We compared these two sets of oscillator strengths under the as- 
sumption of a honaogeneo us abundance, finding that the data 
from Bie mont et al.l (l2Qllh are less consistent internally. Since 
the lines giving discrepant results are not distinguished by their 
strength or excitation potential, this indicates a real limitation of 
the atomic data rather than an impact of possible stratification. 
The list by Biemont et al. also lacks information on two other 
lines, especially the high-excitation Yii 4465 A, which turned 
out to be very useful for stratification analysis. 

On the other hand, all oscillator strengths of Ti lines 
were provided by V ALD and came from a single source 
(Pic kering et al.|[2Q0]]) . The oscillator strength corrections pro- 
vided in the last column of Table \5\ were not applied to Y 
and Ti lines in the stratification analysis. To compute stratifica - 
tion models, we used the IDL code DDAFit ( Kochukhovll20Q7h . 
which provides an interface to the spectral synthesis program 
SYNTH3. DDAFit assumes a step-function parameterisation of 
the stratification profile. In this study we assumed a fixed sharp 
transition -0.2 on the logrsooo scale -between the under- and 
overabundance zones and set the element abundance in the lower 
atmospheric to log(A/^ei/Mot) = -15 dex. The abundance in the 
upper atmosphere was determined by the code by fitting ob- 
served profiles for a predefined set of abundance step positions 
between log T5000 = -2 to -5. Then the fit corresponding to each 
of the stratification models was compared to the one obtained 
with a homogeneous vertical distribution. 

In Fig.[T3]we illustrate the line profiles fits obtained for Y 11. 
The fit provided by an extreme chemically stratified model yields 

^ http://kurucz.harvard.edu/atoms/390 1/ 
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Fig. 8. Surface maps of Y, Sr, Ti, and Cr derived for ip Phe. The star is shown at four equidistant rotational phases with rotational axis 
oriented vertically. Inclination of the rotational axis is 65°. The abundance scale relative to the Sun is plotted on the right. Darker 
regions correspond to a higher relative abundance. The contours for Y plotted with a step of 0.33 dex, 0.13 dex for Sr and Ti, and 
0.04 dex for Cr. 



worse agreement with the observations than does the model with 
a homogeneous Y distribution. 

Standard deviation corresponding to diflTerent stratification 
profiles is reported in Table [6l The first row in this table corre- 
sponds to a homogeneous distribution. There is a marginal indi- 
cation of Y stratification starting at logrsooo ~ -2. This is sug- 
gested by a slightly better fit of the observed line profiles, shown 
in Fig. [HI 

Ti line profiles were examined in the same way as Y. 
Figure [15] illustrates the individual line profile fit in the format 
similar to Fig. [131 The assumption that all Ti is concentrated 
in the upper layers of the stellar atmosphere gives significantly 
worse results, compared to the homogeneous model. Similar to 
the case of Y, we found a marginal indication of Ti stratification 
with a step position at logrsooo = -2. 



7. Summary and discussion 

Our study shows that (p Phe does not possess a strong mag- 
netic field. Owing to the quality of our spectropolarimetric data, 
we were able to set a sensitive upper limit of 4 G on the 
mean longitudinal magnetic field. Our results also do not sup- 



Table 6. Best-fit stratification models and corresponding stan- 
dard deviations for Ti and Y. 



step at 


Ti 






Y 


logTsooO 


^upper 


cr(%) 


^upper 


cr(%) 




-6.55 


0.536 


-7.47 


0.693 


-2 


-5.87 


0.523 


-6.88 


0.648 


-3 


-5.22 


0.680 


-6.22 


0.667 


-4 


-4.42 


0.941 


-5.44 


0.757 


-5 


-3.43 


1.192 


-4.45 


0.874 



Notes. The first row corresponds to a homogenous element distribution. 



port the hypothesis that magnetic field is associated with spots 
of chemical elements. From the spectral lines of inhomoge- 
neously distributed chemical elements, we inferred the upper 
limit of 8- 15 G for {Bj). As in the case of othe r HgMn stars 
studied with high-precision spectropolarimetry ( Auriere et al.l 
2010; Makaganiuk et al. 201 lb), the results of this work clearly 
demonstrate the absence of global magnetic field structures. 
Analysis of LSD profiles also rules out the presence of complex 
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1.2 ~- 




Wavelength 

Fig. 14. Same as Fig. [131 but for the best-fit stratification model 
with a transition step at logrsooo ~ -2. 



magnetic fields, similar to those found in active late-type stars 
(iDonati & Landstreetl l2009h . 

The full rotational phase coverage of our data enabled us to 
confirm the vari ability in the s pectral lines of Y, Sr, and Ti, found 
earlier by Br iquet et al.l (120 10*). These authors also suggest a pos- 
sible variability in Zr. We found a marginal indication of the vari- 
ability in the Zr ii line at 4149.2 A. However, an insuflftcient S /N 
of our data in that spectral region, as well as a small number 
of visible Zr lines, does not allow us to conclusively investigate 
variability for this element. At the same time, a combination of 
superior data quality and the application of LSD technique en- 
abled us to detect a very weak variability of Cr lines. 

We applied the Doppler imaging technique to interpret the 
line profile variability of Y, Sr, Ti, and Cr. The main features of 
the inhomogeneous distribution of all these elements is a com- 
mon spot located on one hemisphere of the star. This spot does 
not occupy the entire stellar hemis phere as in the case of an - 
other HgMn star 66 Eri A studied by lMakaganiuk et all (l2011ah . 



Possibly, the morphological diflTerences of chemical spot distri- 
butions observed in 66 Eri A and (p Phe are related to the absence 
of a close secondary component in the latter star. 

Besides the large primary spots, we found two smaller ones 
for Ti and Cr. These chemical elements show a tendency to con- 
centrate near the stellar equator, unlike Y and Sr, which show lat- 
itudinally extended spots. The spot morphology in HgMn stars 
is poorly studied due to the very small number of stars to which 
DI was applied so far. At the moment, there are only four HgMn 
stars, for which surface maps have been reconstructed (a And, 
AR Aur A, 66 Eri A, and (p Phe). By looking at their maps, one 
can note a tendency for chemical elements that are heavier than 
the mass number Z ~ 37 to have more pronounced spots. These 
chemical elements also show a more contrast between the over- 
and underabundance zones. 

We cannot directly confirm the evolution of chemical spots 
on (p Phe based on co mparing our maps with those published 
by lBriquet et al.l (l2010l) . This is due to the diflTerences in the data 
quality, the time span of observations, and the mapping approach 
used by Briquet et al. In particular, we used multiple spectral 
lines for Doppler imaging, which yields more robust surface dis- 
tributions. Our surface maps are smoother and show less con- 
trast, especially for Ti. In contrast. Briquet et al. finds a number 
of relatively small over- and underabundance spots of Ti. This 
may be related to insuflftcient quality of their spectra or incor- 
rectly chosen regularisation parameter. Briquet et al. did not pro- 
vide any individual line profile fits for Sr and Ti, which does not 
allow us to assess the reasons for the observed diflTerences be- 
tween our maps and theirs. However, we note that there appears 
to be a common overabund ance zone at phases 0.8-0.9 in the 
Ti and Y maps obtained by lBriquet et al.l (l2010l) . The common 
spot in our maps appears at phase 0.2. This diflTerence may come 
from the diflTerently assumed zero phase or insufliiciently precise 
rotational period. 

The evolution of chemical structures claimed by Briquet et 
al. generally consists of a redistribution of small-scale features. 
Our inversion results lack such local features entirely, suggesting 
that the spot evolution in (p Phe may be an artefact of their poorly 
constrained Doppler reconstructions. The primary evidence for 
surface structure variation - the line profile and equivalent width 
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diflTerence at similar rotation phases -has to be examined fo r 
(f Phe in the same way as for a And by'Kochukhov et all (l2007l) . 

It was suggested that heavy chemical elements, such as mer- 
cury, may a ccumulate in th e uppermost atmospheric layers of 
HgMn stars (iMichaud et al.|[T974^) where they could be affected 
by very weak magnetic fields (lAlecian et aDl2Qlll) . We tested 
this hypothesis for Ti and Y with the help of stratification analy- 
sis, assuming that these elements are concentrated only in the up- 
per layers of the atmosphere of ip Phe and are completely absent 
in the lower layers. This extreme stratification model clearly dis- 
agrees with observations. Thus, we conclude that the horizontal 
element inhomogeneities in ip Phe are not associated with high- 
lying chemical clouds. Previous strati fication studies of HgMn 
stars (ISmithI 119961: iThiam et al.ll2QlQh also found no evidence 
of any strong outward abundance increase. Nevertheless, in our 
analysis of (p Phe the formally best fit is achieved for a stratifica- 
tion model with the transition between the over- and underabun- 
dance zones occurring at logrsooo = -2. This stratified distribu- 
tion provides a marginally better description of the observed Ti 
and Y line profiles than a homogeneous abundance. 

The absence of a strong magnetic field that could be asso- 
ciated with the formation of chemical spots suggests the op- 
eration of an alternative process of spot formation. That we 
observe evolving surface inhomogeneities in HgMn stars sug- 
gests a mechanism that aflTects the atomic diflfusion, making it 
time-dependent. Alecian (1998) presents a qualitative scenario 
in which a chemical element is pushed by radiative acceleration 
to the upper layers of stellar atmosphere, forming an unstable 
stratification which changes with time. An increase in concen- 
tration of the element at the lower atmospheric layers absorbs 
photons, thus reducing the radiative acceleration received by the 
same element in the upper layers, leading to cyclic disappear- 
ance and accumulatio n of the upper overabu ndance zone. 

In the recent paper. lAlecian et all (l201ll) has performed a de- 
tailed numerical study of non-equilibrium, time-dependent diflTu- 
sion for a fictitious chemical element similar to mercury. Their 
calculations suggest that the equilibrium s tratification assumed 
by previous atomic diflTusion studies (e.g. jLeBlanc et al.]|2009l) 
may never be attained in reality and that vertical element strati- 
fication may become unstable whenever radiative acceleration is 
equal to gravity. Instabilities first develop in high layers and then 
propagate downwards. The diflTusion time scale can be as short 
as one month at log T5000 = -4, but increases to several years 
in deeper layers. The lack of concentration of Y and Ti in the 
highest layers combined with the claimed spot variation during 
two months seems to be too rapid compar ed to what could be ex - 
pected from the numerical simulations by lAlecian et al. 
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Fig. 9. Comparison of the observed {symbols) and model {solid lines) profiles for Y lines in Phe. The bars in the lower left corner 
indicate the horizontal and vertical scales (0.1 A and 5% of the continuum intensity, respectively). 
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Fig. 11. Same as Fig.|9]but for Ti. 
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